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Abstract 

One of the important goals for future neutrino telescopes is to identify the flavors of astrophys- 
ical neutrinos and therefore determine the flavor ratio. The flavor ratio of astrophysical neutrinos 
observed on the Earth depends on both the initial flavor ratio at the source and flavor transitions 
taking place during propagations of these neutrinos. We propose a model independent parametriza- 
tion for describing the above flavor transitions. A few flavor transition models are employed to 
test our parametrization. The observational test for flavor transition mechanisms through our 
parametrization is discussed. 

PACS numbers: 95.85.Ry, 14.60.Pq, 95.55.Vj 
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Recent developments of neutrino telescopes [l-^] have inspired numerous efforts of study- 



ing neutrino flavor transitions utilizing astrophysical neutrinos as the beam source Q-|25|. 
Given the same neutrino flavor ratio at the source, some flavor transition models predict 
rather different neutrino flavor ratios on the Earth compared to those predicted by the 



standard neutrino oscillations 



(io|. 



In this article, we propose a scheme to parametrize 
flavor transition mechanisms of astrophysical neutrinos propagating from the source to the 
Earth. As will be shown later, such a parametrization is very convenient for classifying 
flavor transition models which can be tested by future neutrino telescopes. 

To test flavor transition mechanisms, it is necessary to measure the flavor ratio of astro- 
physical neutrinos reaching to the Earth. The possibility for such a measurement in IceCube 



has been discussed in Ref. 



,8 



It is demonstrated that the v e fraction can be extracted from 
the measurement of the muon track to shower ratio by assuming flavor independence of the 
neutrino spectrum and the equality of and v r fluxes on the Earth due to the approximate 



v T symmetry 



26 



271 ] . Taking a neutrino source with fluxes of v e and given by 



El e dN u JdE Ve = 0.5E^dN J ,/dE l/lx = 1CT 7 GeV cm~ 2 s~\ which is roughly the order of the 



Waxman-Bahcall bound [281 ]. and thresholds for muon and shower energies taken to be 100 
GeV and 1 TeV, respectively, the v e fraction can be determined to an accuracy of 25% at 
IceCube for 1 yr of data taking, or equivalently to an accuracy of 8% for a decade of data 
taking. Such an accuracy is obtained for a v e fraction in the vicinity of 1/3. The accura- 
cies corresponding to other central values of the v e fraction are also presented in Ref. [8J. 
The to v T event ratio can also be measured in IceCube. However, the accuracy of this 
measurement is limited by the low statistics of v T events. 

Neutrino flavor ratio represented by the ternary plot- To study neutrino flavor tran- 
sitions, we describe the neutrino flavor composition at the source by a normalized flux 



= (0o,e, 0o,At, 4>o,t) satisfying the condition [29 ] 



(fio,a > 0, for a = e,/i,r, (1) 

where each 0o,a is the sum of neutrino and antineutrino fluxes. Any point on or inside the 
triangle shown in Fig. [1] represents a specific flavor ratio characterizing the source. The 
triangular region bounded by vertices (1, 0, 0) T , (0, 1, 0) T , and (0, 0, 1) T contains all 
possible source types in terms of flavor ratios. The pion source and the muon-damped 
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source with flavor compositions $o i7r = (l/3,2/3,0) T and $0,^ = (0, 1,0) T , respectively, are 
explicitly marked on the figure [30]. 

The net effect of flavor transition processes occurring between the source and the Earth 
is represented by the matrix P such that 



o- 



(2) 



where $ 



is the flux of neutrinos reaching to the Earth. We note that our 



convention implies P a p = P{vp —> v a ). 

Q matrix parametrization for flavor transitions of astrophysical neutrinos.- Since the 
triangular region in Fig. [1] represents all possible neutrino flavor composition at the source, 
it is convenient to parametrize $0 by |31 | 

$0 



I Vl + aV 2 + 6V 3 , 



(3) 



where Vi = (1,1, 1) T , V 2 = (0, -1, 1) T , and V 3 = (2,-l,-l) T . Mathematically, Vi/3 
represents the center of the triangle, while aV 2 and 6V3 represent horizontal and vertical 
displacements within the triangle, respectively. The ranges for a and b are —1/3 + 6 < 
a < 1/3 — b and —1/6 < b < 1/3 such that Eq. ([3]) covers all points of the triangular 
region. The pion source and the muon-damped source mentioned in Fig. [1] correspond 
to (a, b) = (—1/3,0) and (a, b) = (—1/2,-1/6) respectively. In general, a source with a 
negligible v T flux corresponds to a = — 1/3 + 6. 

The above parametrization for $0 is also physically motivated. The vector Vi gives the 
normalization for the neutrino flux since the sum of components in Vi/3 is equal to unity, 
while the sum of components in V 2 and that in V 3 are both equal to zero. The vector aV 2 
determines the difference between and v r flux, </>o. M — <po,r, while preserving their sum, 
00,^ + 0o,t- Finally the vector 6V 3 determines the sum of and v T flux, o ,m + 0o,t, while 
preserving their difference 4>o^ — 0o )T . 

Following the same parametrization, we write the neutrino flux reaching to the Earth as 

(4) 



$ = kVi + pV 2 + AV 3 . 



It is easy to show that 



Qn Q12 Q13 
Q21 Q22 Q23 
1 Q31 Q32 Q33 j 



^ /l/3^ 



a 
b 



(5) 
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FIG. 1: The ternary plot for representing the flavor ratio of astrophysical neutrinos. The numbers 
on each side of the triangle denote the flux fraction of a specific flavor of neutrino. The blue 
point, situated on the left side of the triangle, marks the pion source $o,vr = (1/3, 2/3,0) T and 
the red point, situated at the lower-left corner of the triangle, marks the muon-damped source 



0,/i 



(0,1,0)' 



where Q = A PA with 



A 



10 2 
1 -1 -1 

i i -ij 



(6) 



In other words, Q is related to P by a similarity transformation where columns of the 
transformation matrix A correspond to vectors Vi, V 2 , and V 3 , respectively. 
The parameters k, p and A are related to the flux of each neutrino flavor by 

e = K + 2A, M = K - p - A, <fi T = K + p - A, (7) 

with the normalization <fi e + (fi^ + <fi T = 3n. Since we have chosen the normalization o ,e + 
<fiid.ii + <fio,T = 1 for the neutrino flux at the source, the conservation of total neutrino flux 



during propagations corresponds to K — 1/3. In general flavor transition models, K could 
be less than 1/3 as a consequence of (ordinary) neutrino decaying into invisible states or 
oscillating into sterile neutrinos. To continue our discussions, it is helpful to rewrite Eq. (J7J) 

as 

p=((f> T - <f>n) /2, A = e /3 - (0 M + (fi T ) /6. (8) 

It is then clear from Eqs. © and ([SJ) that, for fixed a and b, the first row of matrix Q 
determines the normalization for the total neutrino flux reaching to the Earth, the second 



row of Q determines the breaking of — v T symmetry in the arrival neutrino flux, and the 
third row of Q determines the flux difference (f) e — (0 M + <fi T )/2. 

Compared to P, matrix Q is very convenient for classifying flavor transition models. First 
of all, those models which preserve the total neutrino flux are characterized by the condition 
S a =e pl t P<*P = 1 in the P matrix parametrization. On the other hand, these flux- conserving 
models must give k = 1/3 in the Q matrix parametrization, irrespective of the initial flavor 
composition characterized by parameters a and b. This implies Qu = 1 and Q12 = Q13 = 
from Eq. (151). Clearly the flux-conservation condition in the Q matrix parametrization is 
much simpler. Second, for those models which do not seriously break the — v T symmetry, 
the second and third rows of P are almost identical, i.e., (P Me , P MM , P MT ) ~ (P re , P TA1 , P TT ), 
and the second and third columns of P are also almost identical, i.e., (P eA1 , P^, P T ^) T ~ 
(P eT , Pfir, Ptt) t ■ Using these conditions and the relation Q = A _1 PA, one can show that 
(Q21, Q22, Q23) ~ (0,0,0) and (Q u , Q22, Q32V ~ (0,0, 0) T . Obviously, the approximate 

— v T symmetry is realized by a much simpler condition in the Q matrix parametrization. 
In summary, we have seen that the first and second rows of Q as well as the matrix element 
Q32 are already constrained in a simple way by assuming the conservation of total neutrino 
flux and the validity of approximate — v T symmetry. Hence, under these two assumptions, 
one can simply use the values for Q31 and Q33 to classify flavor transition models. This is the 
most important advantage of Q matrix parametrization. In fact, as will be elaborated later, 
this parametrization is also very useful for discussing the effect of flux nonconservation, i.e., 
the case with k 4^ 1/3. 

n 

It was first discussed in Ref. [7| that the flavor measurement in neutrino telescopes is useful 
for stud y in g neutrino flavor composition at the astrophysical source (for recent studies, see 



Refs. 



31 



32]) and neutrino flavor transition properties during its propagation from the source 
to the Earth (see also Ref. [g]). For probing flavor transition properties of astrophysical 
neutrinos, the authors of Ref. |7|] considered typical astrophysical sources and applied the 
flavor transition matrix P (denoted by \ in the original paper) derived from the standard 
neutrino oscillation model or flavor transition models involving new physics for obtaining 
possible flavor ratios to be measured by terrestrial neutrino telescopes. It was pointed out 
that there are some flavor transition models which can produce rather distinctive neutrino 
flavor ratios on the Earth compared to those produced by the standard neutrino oscillation 
model, even with uncertainties of neutrino mixing parameters taken into account. Hence 
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these flavor transition models can be tested on the basis of their flavor-ratio predictions for 
astrophysical neutrinos arriving on the Earth. In our approach, we test the fundamental 
structure of a given flavor transition model, namely the Q matrix of the model. As noted 
earlier, possible neutrino flavor transition models which conserve the total neutrino flux 
are encompassed by possible values of Q21 and with i = 1, 2, and 3. In the — v T 
symmetry limit, only the values of Q31 and Q33 are relevant. The matrix elements of Q 
can be determined by performing fittings to the flavor-ratio measurements in the neutrino 
telescopes, as will be demonstrated later. The obtained ranges for these matrix elements 
can be used as the basis for testing any flavor transition model. 

Examples.- In the previous section, we have discussed the properties of the Q matrix 
and their advantages. In this section, we shall confirm such properties using a few flavor 
transition models as examples. We begin by considering the standard three-flavor neutrino 
oscillations. It is well known that 

3 

^El^l 2 !^! 2 ' ( 9 ) 

i=l 

for astrophysical neutrinos traversing a vast distance where U a i and are elements of 
the neutrino mixing matrix. It is easily seen that P°^ c = P^ . Because of the probability 
conservation, the flux of neutrinos on the Earth also satisfies the normalization condition 
given by Eq. ([T]). We first compute Q osc in the tribimaximal limit 33] of neutrino mixing 
angles , i.e., sin 2 #23 = 1/2, sin 2 6*12 = 1/3 and sin 2 #13 = 0. In this limit, the — v T 
symmetry is exact, to fact, an exact „„-„, symmetry amounts to the condition \UJ = \U H \ 
for i = 1,2,3 [26|, |27|. This condition can be realized by having both sin 2 623 = 1/2 and 
sin #13 cos 5 = 0, which are respected by the above tribimaximal limit of neutrino mixing 
angles. Denoting P osc in this limit as Pq sc , we have 



f 5/9 2/9 2/9 ^ 



pose 



(10) 



2/9 7/18 7/18 
^2/9 7/18 7/18 J 

Since — v T symmetry is exact in this case, one can see that the second and the third rows 
of Pq sc are identical, so are the second and third columns of Pq sc . The corresponding Q 
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matrix in this limit is given by 



^10 0^ 



Q° sc = A- 1 PS" A 



• (11) 
yo 1/3 J 

As expected, Qq s u = 1 and Qq s i 2 = Qo s i3 = 0- Furthermore, any element in either the 
second row or the second column of Qjf vanishes. 

We can compute the correction to Qq sc as neutrino mixing parameters deviate from the 
tribimaximal limit. We consider such deviations for # 13 and # 23 while keeping sin 2 # 12 = 1/3. 
In this case P osc = Pq sc + P° sc + • • • where P° sc is the leading order correction in powers of 
cos 26*23 and sin 6*13. One has [34 j 

/ n . A 



pa 



e -e 
e -e 
-e e 



(12) 



where e = 2cos2^ 2 3/9 + \/2 sin 6*i 3 cos 5/9 with 5 the CP phase. Taking into account P° sc , 
we obtain Q osc = Q° sc + Q° sc with 



/ 



gosc = A -lposc A 


















-3e 





— e 






(13) 



Therefore Q osc is given by 



A -lposc A 



1 






-3e 



(14) 



y0 -e 1/3 J 

Because of the correction term Q° sc , one can see from Eq. tHM that the v u 



v T symmetry is 

broken since Q 1 ^ and Q^ S 2 C are nonvanishing. Focusing on the third row of Q osc , we obtain 
A = b/3 - ae from Eqs. f|5l) and ffMl). 

We next consider models of neutrino decays. Flavor transitions of astrophysical neutrinos 
due to effects of neutrino decays were discussed in Ref. [6]. The simplest case of neutrino 
decays is that both the heaviest and the middle mass eigenstates decay to the lightest mass 
eigenstate. We first assume the branching ratios for the above two decays are both 100%. 



Elements of subleading matrices Q' dcc and Q'[ 



y/dcc 





12 


21 


23 


32 


Qf ec 


-2(l-r-s)(ei + e 2 )/3 


-(l + r)ei - (l + s)e 2 


rei - (1 - s)e 2 


Hex + e 2 ) - e 2 ]/3 


Q//dec 


2(1 - r- s)ei/3 


(1 + s)ei - (r - s)e 2 


-ei - 2e 2 


-[(l + r-s)ei + 2e 2 ]/3 



TABLE I: Nonzero elements for subleading matrices Q' dec and Q'[ dcc . The indices 12,21,23, and 
32 in the heading of the table denote the positions of matrix elements, r and s denote branching 
ratios for the decays v% — > z^ 2 and v% — > fi, respectively, in the case of normal mass hierarchy, and 
branching ratios for the decays f 2 — > v\ and vi —> z^, respectively, in the case of inverted mass 
hierarchy, ei = (cos 2# 2 3 — \[2 sin #13 cos <5)/3 and e 2 = cos 2# 2 3/2 — t\. To obtain these expressions, 
we have taken sin 2 #i 2 = 1/3. 



Under this condition, the transition matrix is given by P dec 



a/3 



\U a i\ 2 for the normal mass 



hierarchy and P d ^ c 
then reads 



I U a 3 1 2 for the inverted mass hierarchy. The corresponding matrix Q 



dec 



/ 



Q 



dec 



V 



\u -i 2 



1 

3(|f/ Ti | 2 -|cy 2 )/2 

2 1 ITT .I2N 



(15) 



(l^l 2 + l^| 2 )/2 y 

where j '• = 1 for the normal mass hierarchy and j = 3 for the inverted mass hierarchy. One 
can see that Q d l c = 1 and Qf^ = Qf™ = 0. Furthermore, in the limit of exact — v T 



symmetry, one has \U. 



\U, 



i-'j 1 



such that the elements in both the second row and the 
second column of Q dec vanish. If branching ratios for the above decays are not 100%, the 
resulting Q dec matrix would be different but nevertheless gives rise to the same neutrino 
flavor ratio on the Earth. It is interesting to see that all the nonvanishing elements of Q dcc 
are located in the first column. Hence, following Eq. ([5]), the neutrino flavor ratio on the 
Earth is independent of the neutrino flavor ratio at the source in this scenario. 

Let us consider another neutrino decay scenario where only the heaviest neutrino decays. 
Following earlier treatments, we set sin 2 612 = 1/3 while allowing 623 and #13 to deviate from 
7r/4 and 0, respectively. For the normal mass hierarchy, we write Q' dec = Qo dec + Q'i° c where 
Qo dcc is the leading term obtained in the limit sin 2 ^ = 1/2 and sin# 13 = 0, while Q[ dec is 
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the first-order correction which is linear in cos2# 23 and sin 6 is . We find 



Qo 



/dec 



4 + 2(r + s) 2 

1 + s 



- 2(r + s) 


1 - s 



(16) 



and 



/ 



/dec 







\ 



21 



J 23 



V 



o (Qf c ) 32 < ! 



(17) 



/ 



(gf c )i2 o 

(Q[ de ^ 

)32 

where r and s are the branching ratios for the decay modes u 3 — > v 2 and u 3 — >■ respectively. 
The nonzero elements of Q' dec are given in Table I. 

If us exclusively decays into either u 2 or i>i, one has r + s = 1. In this limit, Qffi* = 1, 
Q'f 2 ec = Q'f 3 cc = as expected. One also observes that the elements in the second row and 
the second column of the leading matrix Q'q oc vanish due to z/ M — u T symmetry. Finally, the 
third row of Q' dcc gives rise to A = [(1 + 36) + (1 - 3b)s] /18 + a [s(ei + e 2 ) - e 2 ] /3. Focusing 
on the leading order contributions, one has A > since b < 1/3; i.e., e is either equal or 
larger than (0 M + <p T )/2 irrespective of the flavor ratio at the source. For comparison, the 
standard oscillation scenario gives A = b/3 at the leading order, which is either positive or 
negative depending on the sign of b. 

For the inverted mass hierarchy, we denote r and s as branching ratios for the decay 
modes u 2 — > vi and v 2 — > z/ 3 , respectively. We obtain Q" dcc = <5o' dcc + Q'[ dec with 



t 4 + 2(r + s) 0^ 



Qo 



//dec 



V 



o 

r — s 




2 



(18) 



and 



/ 



Q" 



(Qf ec )i2 o 



\ 



V 



(19) 



/ 





//dec_ (g//dcc)2i Q (g/ , dcc)23 

o (gf cc ) 32 o 

The nonzero matrix elements of Q'{ dcc are given in Table I. In the limit r + s = 1, Q'lf 00 = 1, 
Q" dec = Q'( dcc = as expected. It is also observed that the elements in the second row and 
the second column of the leading matrix <5o' dcc vanish. Finally, the third row of Q" dcc gives 
rise to A = (r - s + 66)/18 - a[(l + r - s)ei + 2e 2 ]/3. 
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As the last example, we discuss neutrino flavor transitions affected by the decoherence 



effect from the Planck-scale physics 
that 



36 



38 



35|. In a three- flavor framework, it has been shown 



pdc 



+ 



+ -e 
6 



-78 d 



(Ul l + Ul 2 -2U%)(Ul l + U 2 a2 



2^ 3 )] 



(20) 



where 73 and 7§ are eigenvalues of the decoherence matrix, and d is the neutrino propagating 
distance from the source. The CP phase in the neutrino mixing matrix U has been set to 
zero. Taking 73 = 7s = 7, we obtain Q dc = A _1 P dc A = Q dc + Q dc where 



Qo 



dc 



/ 



V 



and 



10 ^ 

J 

( 



(21) 



Qi 



dc 



\ 



V 



/ 



00 

-3e , (22) 
-e 

with e = 2 cos 2^23/9 + \/2 sin 6*13/9. From the definition right below Eq. ( fT2|) . we note that 
e = e(5 = 0). One can see that Qn — 1 ; an d Q12 = Q13 = 0. Furthermore the elements in 
the second row and the second column of the leading matrix Q^ c vanish. In the absence of 
the decoherence effect, i.e., 7 — > 0, it is seen that Q dc reduces to Q osc . In the full decoherence 
case, i.e., e~ ld — > 0, we have k=1/3, p = X = such that e : M : <p T = 1 : 1 : 1. 

Probing Q by measuring flavor ratios of astrophysical neutrinos - We have shown that 
the flavor transitions of astrophysical neutrinos can be parametrized by the matrix Q. As we 
have argued earlier, the Q matrix is very convenient for classifying flavor transition models. 
One could determine the matrix elements Qij by measuring flavor ratios of astrophysical 
neutrinos arriving on the Earth. In this regard, we derive from Eqs. ([5]) and ([8]) that 

1 

-( 

3 

1 

-( 

3 



3(/ T (a,6)-/ A ,(a,6))/2 
/ e (a,6)-(/ M (a,6) + / T (a,6))/2 



-Q21 + aQ 22 + bQ 2 3j /«(a,6), 
-Q31 + aQ 3 2 + bQ 3 3) /«(a,6), 



(23) 
(24) 
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where f a is the fraction of u a , i.e., f a = </> a / (0 e +</v+0i-) = q /3k. In the above equations, we 
have explicitly denoted the dependence of f a on the source parameters a and b. Furthermore 
we also indicated that k is generally a function of source parameters since the total neutrino 
flux is not necessarily conserved during neutrino propagations. 

In the flux-conservation case, Qu = 1 and Q12 — Q13 — 0, which gives k = 1/3. In 
principle, the matrix elements Qu in the second row of Q can be solved from Eq. (1231) by 
inputting three sets of f a (a, b) measured from three different astrophysical sources. Here we 
assume precise knowledge of parameters a and b from each source. The matrix elements Q 3 i 
in the third row of Q can be solved from Eq. (j2"31) in a similar way. In the case that — v T 
symmetry is not significantly broken, one expects Q211Q221 and Q23 are all suppressed. 
Therefore, it is more involved to probe the second row of Q than to probe the third one. To 
probe the latter, we have 

f e (a, 6)/3 - (/>, b) + f T (a, b)) /6 « ±Q 31 + bQ 33 , (25) 

since Q32 is also suppressed due to the approximate — v T symmetry. We note that f a (a, b) 
on the left-hand side of Eq. fl25|) only depends on b. It is possible to solve for Q 3 i and 
Q 33 if the measurement on f e — (/ M + f T )/2 can be performed with respect to two different 
astrophysical sources where the value of the b parameter in each source is known. 

In the case of flux nonconservation, the function n(a, b) is not known since it is difficult 
to determine the absolute flux of astrophysical neutrinos at the source. Hence one cannot 
directly solve for Q 2i and from Eqs. (123]) and by inputting f a (a,b) from measure- 
ments. On the other hand, the signature for k 7^ 1/3 could still be detected by the following 
consistency analysis. We recall from Eq. (12^1) that the third row of Q is related to the 
measurement by 

f e (a, b) - (f^a, b) + f T {a, b)) /2 « Qq 3 i + bQ 33 ) /«(o, b). (26) 

As it was just argued, one could set k = 1/3 in the above equation and invoke two astro- 
physical sources to solve for Q 3 i and Q 33 . However, taking this set of Q31 and Q 33 as an 
input, one expects that the right-hand side of Eq. (T2"6"|) is likely to be inconsistent with the 
left-hand side obtained from the third astrophysical source. 

It is clear that the knowledge of the neutrino flavor ratio at the source is crucial for probing 



the matrix Q. Previous studies 



39|, 



40| pointed out that this ratio is energy dependent for 
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a general astrophysical source. For parent pions with an E~ 2 energy spectrum, the flavor 
ratio of neutrinos arising from the decays of these pions and the subsequent muon decays is 
0o,e : 00,^ : 0o,r = 1 : 1.86 : at low energies [411] where energy losses of pions and muons 
in the source are negligible. The ratio 0o,e/</>o,/i however decreases with the increase of 
muon (pion) energy and eventually approaches zero. This behavior results from the above- 
mentioned energy losses which are important at higher energies. Recently, a systematic 
study on possible neutrino flavor ratios from cosmic accelerators listed on the Hillas plot 
was initiated [42]. The neutrino flavor ratio at the source depends on the spectrum index 
of injecting protons, the size of the acceleration region, and the magnetic field strength at 
the source. In some regions of the above-mentioned parameters, the neutrino flavor ratios 
are energy dependent, while in some other parameter regions they could behave as those 
of a pion source or those of a muon-damped source, which are both energy independent. 
In the following, we illustrate the determination of Q31 and Q33 by measuring flavor ratios 
of astrophysical neutrinos arriving on the Earth from a pion source and a muon-damped 
source, respectively. 

To determine Q31 and Q33, we assume an exact — v T symmetry so that M = <f> T . 
The measurement of muon track to shower ratio §J in a neutrino telescope such as IceCube 
can be used to extract the flux ratio R = M /(0 e + <f> T ). Clearly R depends on the source 
parameter b and the matrix elements Q31 and Q33 as can be seen from Eq. fl25|) . One can in 
principle disentangle Q31 and Q33 by measuring R from two sources with different b values, 
say a pion source with 6 = and a muon-damped source with b = —1/6. Taking into account 
experimental errors in determining R, the ranges for Q31 and Q33 in a given confidence level 
can be determined by the formula 

2 ( ^7r,th Rjr,exp \ ( Rp,th -^/i,exp 



^2 I - ll, 7r,tn -"TTjexp 1 _|_ < -^"fjixa ^"/i.cxp 1 (27) 

,exp / 

where R Wjex p and -R^exp are experimentally measured flux ratios for neutrinos coming from 
a pion source and muon-damped source, respectively, while R^^h and R^th, which depend 
on Q 3 i and Q 33 , are theoretically predicted values for R T and R^ respectively. Furthermore, 
0Xr,e X p = (A#7r/#7r)K,exp and a Rfi exp = ( AR^ / R^) R^ >exp with AR n and AR^ the experi- 
mental errors in determining R for neutrinos coming from a pion source and muon-damped 
source, respectively. One does not need to include uncertainties of neutrino mixing angles 
9ij and CP phase 5 in Eq. ff27l) since their effects are already embedded in Q 31 and Q 33 . 

12 




FIG. 2: The fitted la (solid line) and 3a (dashed line) ranges for Q31 and Q33. The left panel 
is obtained with measurement accuracies AR n /R n = AR^/R^ = 10%, while the right panel is 
obtained with AR^/R^ = AR^/R^ = 20%. The circle describes the best-fit parameter values 
Q31 = and Q33 = 0.33, corresponding to the input flavor transition model. For reference, the 
parameter values for the neutrino decay scenario given by Eq. f)15j) are denoted by the triangle 
and the square, respectively, for normal and inverted mass hierarchies. The former corresponds 
to (Q3i,Q33) = (0.5,0), while the latter corresponds to (Q3i,Q33) = (—0.5,0) for neutrino mixing 
parameters taking the tribimaximal values. 

Let us first take the input (true) flavor transition mechanism to be a standard neutrino 
oscillation model with neutrino mixing parameters taking the tribimaximal values. One 
expects that -R^exp and R^ exp are around 0.50 and 0.64, respectively. Applying the x 2 
analysis, Eq. fl27|) . with given accuracies (Tr^ and <tr p , the fitted la and 3a ranges for 
Q31 and Q33 are presented in Fig. [2j We note that the left panel is obtained with AR n / R^ = 
AR^/R^ = 10% while the right panel is the result of taking AR n /R 7T = AR^/R^ = 20%. 
For both measurement accuracies, the neutrino decay scenario given by Eq. (1151) can be 
ruled out at the 3a level for both mass hierarchies. We stress that the confidence ranges in 
Fig. [2] can be used to test any model with specific values for Q31 and Q33. 

We next consider the case where the input flavor transition model is the neutrino de- 
cay scenario given by Eq. ( jl5j) with normal mass hierarchy. This model corresponds to 
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FIG. 3: The fitted lcr (solid line) and 3a (dashed line) ranges for Q31 and Q33. The left panel 
is obtained with measurement accuracies AR n /R n = AR^/R^ = 10%, while the right panel is 
obtained with AR^/R^ = AR^/R^ = 20%. The triangle describes the best-fit parameter values, 
(Q31JQ33) = (0.5,0), corresponding to the input flavor transition model. The circle corresponds 
to the standard neutrino oscillation model, while the square corresponds to the neutrino decay 
scenario given by Eq. (|15p with inverted mass hierarchy. 

{Q311Q33) = (0.5,0) for neutrino mixing parameters taking the tribimaximal values. Hence 
one expects that R n ,ex P and R^ exp are both around 0.2. Applying the x 2 analysis, we ob- 
tain the fitted la and 3a ranges for Q31 and Q33 as shown in Fig. |3j Once more, the 
left panel is obtained with AR^/R^ = AR^/R^ = 10%, while the right panel results from 
ARx/Rn = AR^/R^ = 20%. For both cases, it is seen that the standard neutrino oscillation 
model and the neutrino decay model with ((531,^33) = (—0.5,0) (inverted mass hierarchy) 
can be ruled out at the 3a level. 

Finally, if the input flavor transition model is the neutrino decay scenario given by Eq. ( fl5l) 
with inverted mass hierarchy, i.e., (Q3i,Q33) = (—0.5,0), one expects that R W)Cxp and -R^exp 
are both around 1.0. Applying the \ 2 analysis, we obtain the fitted la and 3a ranges for 
Q31 and Q33 as shown in Fig. HJ For AR n / R v = AR^/R^ = 10% (left panel), it is seen that 
the other two models displayed on the figure can be ruled out at the 3a level. However, for 
AR^/R^ = AR^/R^ = 20% (right panel), the standard neutrino oscillation model cannot 
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FIG. 4: The fitted la (solid line) and 3a (dashed line) ranges for Q31 and Q33. The left panel 
is obtained with measurement accuracies AR n /R n = AR^/R^ = 10%, while the right panel is 
obtained with AR^/R^ = AR^/R^ = 20%. The square describes the best-fit parameter values, 
((531,(533) = (—0.5,0), corresponding to the input flavor transition model. The circle corresponds 
to the standard neutrino oscillation model, while the triangle corresponds to the neutrino decay 
scenario given by Eq. (|15p with normal mass hierarchy. 

be ruled out at the same confidence level. 

Conclusion. -In summary, we have proposed to parametrize the flavor transitions of propa- 
gating astrophysical neutrinos by the matrix Q, which is related to the usual flavor transition 
matrix P by Q = A~ l PA where A is given by Eq. ([6]). We have argued that it is much 
easier to classify flavor transition models by the Q matrix parametrization, where each row 
of Q carries a clear physical meaning as illustrated by Eq. ([5]) . We have also argued that the 
signature for flux nonconservation might be detectable if it is possible to observe sufficient 
numbers of astrophysical neutrino sources with different flavor ratios. For the case of flux 
conservation, the above observations can probe the second and the third rows of matrix Q 
in a model independent fashion. 

For illustration, we considered the determination of the Q matrix in the exact — v T 
symmetry limit. The relevant matrix elements in this case are Q31 and Q33. We proposed 
to determine them by measuring the flux ratio R = <p^/ (<f) e + <ft T ) for astrophysical neutrinos 
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coming from a pion source and those coming from a muon-damped source respectively. 
We fitted Q31 and Q33 to the measured flux ratios -R^exp and R^ exp using Eq. (127]). The 
ranges for Q 31 and Q 33 are presented up to the 3a confidence level for three different input 
models for neutrino flavor transitions. We have found that the measurement accuracies 
AR^/R^ = AR^/R^ = 10% are sufficient to discriminate among the standard neutrino 
oscillation model and neutrino decay scenario given by Eq. (|T5|) for normal and inverted 
mass hierarchies. We reiterate that the confidence ranges in Figs. El [3] and H] can be used to 
test any flavor transition model with specific values for Q31 and Q33. 

Taking a neutrino source flux E 2 U dN u JdE Ue = 0.5£ 2 dK.JdE^ = 10" 7 GeV cm" 2 s _1 , 



which is roughly the order of the Waxman-Bahcall bound [28], the accuracy AR/R = 10% 
is reachable by a decade of data taking in Icecube fl, as stated in the be ginning of this 
article. However, we stress that the Waxman-Bahcall bound is for diffuse neutrino flux. 
The flux from an individual point source is smaller. Hence it could take more than a decade 
to reach a 10% accurate measurement on R arising from a point source. The radio extension 
of IceCube jsj is expected to accumulate neutrino events at a much faster pace. It is crucial 
to study the efficiency of flavor identification in this type of detector. 

In this work, the Q matrix is probed by assuming an exact — v T symmetry and a precise 
knowledge of the neutrino flavor ratio at the source. Away from the — v T symmetry limit, 
the second row of Q and Q32 shall become relevant in addition to Q31 and Q33. Furthermore, 
the statistical analysis outlined by Eq. ([27]) should be refined once the uncertainty of the 
neutrino flavor ratio at the source is taken into account. We shall address these issues in a 
future publication. 
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